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A THEORETICAL INVESTIGATION OF LONGITUDINAL STABILITY OF AIRPLANES WITH FREE 
CONTROLS INCLUDING EFFECT OF FRICTION IN CONTROL SYSTEM 

By Hauky Ghkknberg atul Lkonard Sternfikld 



SUMMARY 

The f'f'/atlon hvtweeii the drmtor huKje-mometd paranieUrs 

(I fid the CO /it rot forces for changes in forward speed and In. nianeu- 
vers is shown for sereral values of static stabiHty and elerator 
mass balance. 

The stability oj the short-period oscillations is shown, as a 
series of boundanea tjinng the limits oj the stable region in 
terms of the elevator hinge-moment parameters. The ejfects 
oj static stability^ elevator moment oj inertia^ elevator mass 
unbalance, and airplane density are also considered . Dynamic 
instability is likely to occur ij there is mass unbalance oj the 
elevator control system combined with a small restoring ten- 
dency {high aerodynamic balance). This instability can be 
prevented by a rearrangement oj the unbalancing weights which, 
however, involves an increase oj the amount oj weight necessary, 
ft can also be prevented by the addition oj viscous jriction to 
the elevator control system lyrovided the airplane center oj 
gramty is not behind a certain critical position. 

For high values oj the density parameter, which correspond 
to high altitudes oj Jlight, the addition oj moderate amounts of 
viscous jriction may be destabilizing even when the a irplane is 
statically stable. In this case^ increasing the nscous jriction 
makes the oscillation stable again. The condition in which 
viscous friction causes dynamic instability oj a statically stable 
airplane is lim ited to a definite range oj hinge-moment param- 
eters. It is shown that, when viscous jriction causes increas- 
ing oscillations, solid friction will produce steady oscillations 
having an amplitude proportional to the amount oj jriction. 

INTRODUCTION 

Tho clffu'ts of aerodynamic Imlttiun* and mass unbalance 
of th(^ t'lovjitor on the dynamic stal)i]ity of tlie airplaiu^ ai*e 
(liscusscd in a previous report on cotitrol-free stability (r(»fer- 
eiice 1). It was foinid tlieoretieally in ref(M-enee I and ViM'ilied 
ill iliglit (^refereiiee 2) tluit, if the elevator is too closely 
balance*! aerodyuamicaliy and has a siiflicient amount of 
mass uiil)ahuiee (which tiMids to depress the elevator), in- 
ereasini^ oscillations of short ])(»riod may occur. Other lliii;lit 
tests (reference i)) showed, however, that mass unbalance of 
the elevator control system improves the static stability of 
an airplane, that is, increases the slope of the curve of stick 
force against speed in level flight and of the curve of stick 



j force against normal acc(»leration in maneuvers, Subse- 
i (pient work (reference 4) has indicated that a control surface 
I with ])ositive lloatitifr tendency (tendency to float ai^ainst 
i the relative wind), when used as a rudder, is ell'ective in im- 
I provinir control-free static stability. A theoretical analysis 
j (reference o) showed that a rudder havini; a ])osilive lloatinj: 
ratio may, under the influen<*e of solid friction in the control 
system, build up steady oscillations of the airplane and 
rn<lder. These steady oscillations have been observed in 
(lij^ht t(»sts (reference ())- These results suggested an investi- 
gation of the behavior of an airplane e{juipj)ed with an 
elevator having a positive floating tendency. This type of 
elevator was not considered in any of the previous investi- 
gations. 

The ])urpose of th(» present re])ort is to make a theoretical 
analysis of the control-free longitudinal stability of an air- 
plane, which takes account of this current trend toward a 
posit iv(i floating t(»ndency in control-surface design and 
covers, in general, a mtich wider range of ])aramcters than 
the investigation of r(*ference I. These i)arameters include, 
for the el(»vator control system, restoring tendency, floating 
tendency, mass unbalance (hobwTight control), moment of 
inertia, and viscous and solid friction and, for the airplane, 
ihMisity and cetiter-of-gravity position. 

The method of analysis of dynamic stability is based on 
the classical tluM)ry of Bryan and Jiairslow extended to iu- 
I elude movements of the controls and tluMr couplings with 
i the airplane motions. Friction is treated in the same way 
as in the approximate m(»tho<l of reference 5, in which solid 
friction is replaced by an eciuivalent viscous friction. 
I Before the analysis of dynamic stability is presented, some 
j discussion is given of the ellect of the various parameters on 
the elevator forces for trim and for acceleration — charac- 
i t eristics considered important to Hying (pialtitics. The 
; stability of the short-period oscillations, with and without 
I friction in the control system, is then consi<lcred. The effects 
of wj'ights addt^l to the system to modify the static and dy- 
namic stability are discussed. The trends to be expected 
are illustrated by a series of calculations and charts based 
on a typical airplane. The stability of the long-period 
(phugoid) oscillations is not discussed because of its relative 
unimportance. 
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SYMBOLS 

wiiiiT nspi'Cl iMlio 
mil aspect vtxuo 

cocdiciiMits ill siahiliiy iM|UjUi(Mi 
winir span 

olcvalor liiniro-tnomtMU ('(x llicirni 



fricl iomil 



{;J's..) 



(•(M'Hicicrii 



iif 



411. 



A,.= 



pS,c,r 
SA__ 

/•r 



-V, 



applied liiiifriMHonuMil cocirH itMit 
iiirplant! lift (MH-flincMU ( 
lift coolliciont of tail 

pitchinjr-rnonuMit cot^flicitMU ahoiil aiipljiiu' 

center of srravity 
win?; chord 
elevator chord 

difFeroiitial operator ( ^^-^ 

constant term in stability ccpiation 
stick force; imsitivc lor pull 

/(I FA 

stick- force gradient in nianeuv<M-s ^^^^ J 

stick-force j^radient for level tlii^dit 
acceleration of gravity 

hinge moment; positive when tends to de- 
press trail ing edge 

mass moment of elevator ahont its hinge; 
positive when tailheavy 

mass moment of control stick about its ])ivol ; 
])Ositive when stick tends to move forward! 

frietional hinge moment 



N 



moment of inertia of elevator about its 
hinge 

moment of inertia of <*ontrol stick aboiit its 
pivot 



radius of gyration of airplane about F-axis 



distance between airplatie center of gravity 
and elevator liinge 



AT 



V 

w 

X 



leiiglh of control si irk 

pitching moment about air|)hine ciMiler ol 

gravity 
mass of aii plaiK^ 

number of cycles re(|un*e(i foi- oscillation to 
damp to half amplitude 

normal a<'ce[(M'alion \)cv (j of airplane due to 
curvature of ilight path; aceeleromeler read- 
ing minus com])onenl of gravity lorce 

period of oscillation, seconds 

tlynamic pressure 

ehn^ator area 

tail area 

wing area 

distance in half-chords VlVtjv) 

time reriiiired for os<Milation to <lamp to ludf 

amplitmle, s(»conds 
time 



forward velocity 

change in forward velocity from trimmed value 
weight of airplane 
longitudinal forc(-. jiositive forward 
distance of center of gravity from aero<lyna- 
mic center; positive when center of gravity 
is ahead of aerodynamic center 
normal force; positive downward 
angle of attack 
angle of attack at tail 

deflection of elevator; positive for downward 

motion of trailitig edge 
amplitude of elevator oscillation 

angle of down wash 
control gearing {BJhA 
angle of pitch of airplane 

deflection of control stick; positive for for- 
ward motion of stick 
c(miplex root of stability eijuation 
f real and imaginary i)arts, respectively, of \ 

airplane-(i(Misity parameter imipSJi) 
mass densitv of air 
\Vhenev(U- \\ a, a, o] 6, IM. DO, Db, and IPa are used 
as sid)scrii)ts, a derivative is indicated. For example. 
V -^^^ antl C, Whenever a dot is used above a 



dm 



1.= 



symbol, it denotes diflerentiation with respect to time. 
All angles are mt^asurcd in radians. 

METHOD OF ANALYSIS 

Four degrees of freedom— forward speed, angle of attack, 
angle of l)itch, and elevator denection-are generally in- 
volved in tlu^ problem of (;ontrol-fr(»e stability. To each 
clegree of freedom, there corresponds an e(iuation of equihb- 
rium Ix^tweeii inertinl ami aerodymimic forces or moments. 
By us(^ of wind axes, the fotu- e(|uations become, for level 
Ilight, 



rifKoiurrrrAL i wksiicatiox ok L(»\(;rrri)i \al ^lAitiLiTV or aikplaxks with ihkk ('o\Tia)Ls 



IIy\ V • !Lct -t- ///.« //..« - mo I- /Ao. 

wliich vww 1)0 written in iioiuiinuMisiotuil Umw ns 



(1; 



la applying ecjuations (i) to dynjiniic staliility, certain approximations may Ix' madt'. For inslanci*, sliorl-ptM'iod oscil- 
lations (of tlio order of 1 sec) involve ne.iriij?ibl(» ciianires in forward speed, whi<'ii nuiy therefore l)e neglected in studying the 
sliort-period oscillations. In fact, the |)(»riod and dam|)ing of these oscillations cnn he ohtained to a high degree of accuracy 
by using only the last three of the eciuations (1) and setting u = 

Equations (1) then become 



By setting 



+ -r{(\,,-'>Arr(JiLYl))IJd i-iC,„,~\-C\^,D)5, =i) 



(2) 



it can l>e shown (refererwe 7) that X must Ik* a root of a (pnirtic efpuUion formed hy writing 



'2A„fll^y-\ 



= 0 



The resulting stal)ility ecpiation may h<^ w ritten as 

where .1, B. (\ I'J, and /' arc* functions of the stability d(»rivativ<'s. 



(3) 



The stutly of tlu» effects of (hirerent parameters oti the 
control-free stability was made by a s(»ries of co!nputati(ms 
for an average airplane* having tiie charactcnslics given 
hereinafter. The current trend toward a positive floating 
tendency in control-surface design suggests the use of C,,^^ and 
T/,, as th(* fundamental variables to Ix* used in expressing 
stability and control characteristics. The results are pre- 
sent(Ml as a sei*i(»s of figur(*s that show the Halations belwiM'n 
(\ and (\, which, with the other derivatives lixed. satisfv 
the conditions for neutral dynamic stability aiid neutral 
static stal)ility. 

A curve for neutral dynamic stal)ility is the l>oundarv 
divichng the region of increasing oscillations from the region 



of damped oscillations and is obtained from Koutlfs 
discrimimmt 

/^(?A;-.1/l-- 
Th<* corulition for neutral static stability is that 

F = {) 

Tlui stability e(i nation (:]) has four roots. A pair of 
com[)h»x roots indicates au oscillatory mode and a real root 
iiulicates an aperiodic mode. The real part of the complex 
I'oot tlelermines the damping; the imaginary part determiners 
the period of the oscillations. More spocilicaliy. if there is 
a |)air of complex roots 
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iho poriod in sccotuls is irivon l)y 



21' 



and the Uino in socojids to damp to lialf ampliliidc isi^ivtMi l>v 

'/',«=- 



For Hii airplano at ('onstant speed, their may l)e two 
oseillatoi'v modes, Ihovi' may be only one oseillatory mode, or 
tiic motion may ho entirely aperiodic. In cases i[i winch 
there nvc oscillatorv' components, one of the oscillations may 
l)e poorly damped and (?ven l»eeome unstal)io. 

Tlic averairo nirplnno (ni which the calculations of this 
report are hnsed is of conventional desifrn. The char- 
acteristics of the airplane arc 

- - . <) 

k, . 1. 5 

^"^J^,. . 3.3 



I., ft- 



] 

_ 0. IS 

2 

- - - 0. 55 

-_ I. 5 

Tlic basic stability derivatives and parameters obtained from 
these airplane characteristics l)y methods shown in appendix A 
are 

(\-- 

'\ 

0. 486 



('™„. -<' 07 



(witli ao friction) 
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3. 8 
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-8. 9 


0. 486 


^v----- 


3. 22Ck 


-1. 54 





23. 2 


-0. 97 


^ — 


-I0..'-,5C.„ 


-I 







The followinji parameters of the airplane were varietl: 

fi airplane-density parameter 

C,„^ control-fixed static-stability parameter 

The followini^ parameters of the elevator control system 
were varied: 

Ch^^ floatinji: tendency 

Caj restoring tendency 

Cf,^^ elevator-dampinji; parameter 

moment-of-inertia parameter of elevator about its 
hinge 

i, moment-of-inertia parameter of control stic^k about 
its liin^c 

h mass-moment parameter of elevator control system 

about elevator hinge 
/tp mass-moment parameter of elevator alone about its 

liinge 

As has been pointed out, the stability boundaries were 
plotte<l, in most cases, in terms of C^^^ and as the varia- 
bles of tiic coordinate system. In analyzing the eflects of 
friction in the control system, C^^ and Chj,^ were used as the 
plotted variables in some figures whereas C^^^ and C,,^ were 
used in others. The effect of the otlier parameters is found 



by varying them one at a time, tln*ough a range of values, and 
showing for each parameter a s(M-ies of stability boimdaries. 
The size of the airplane, wing loadiuir, and altitude are 



combined in the parameter ^. which is 



A variation 



in M thus could be tluo to a variation in size, wing loading, 
or altitude, or any combination of these. The range of values 
of M cov(Ted in t he prescMit report and some typical correspond- 
ing values of wii^g loading, altitude, and size are given in the 
followin<: table: 
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Sc:i lovt'l 
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Merlin 
\\ ini; 

Ul.i 
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The range of C,„^ and the corresponding centcr-of-gravity 
positions are as follows: 





1 

(fmction nf ; 
iiintn wing | 


-0. 232 
0 

. -'32 


(1.05 [ 
-.05 ' 



The ranges of values of the other parameters, for a small 
airplane (chord, 7 ft), are as follows: 
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STATIC STABILITY AND RELATION TO CONTROL FORCES 

The connection between the static stability and the 
airplane and control parameters is established to assist 
in the ititerpretation of the results obtaiiunl hereinafter. 
Inequations (1) can be applied to static stj,ibilitv by setting 
all terms containing D and D- equal to zero and solving the 
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r esult equal ions sitnultaiu'ously for the vjiriation iti 
lonviinl s|)('(Hi witii an appli(Ml rlovator iiini^o nioiiKMit. 
For iovrl lliirht. 0 is also zcm-o and the ivsultinir o(inations tivv. 

^'r^t i- .> ' =:{) 

Till* variation of sticlv forrr with fractional cluuifre in forward 
speed is 



As shown m a[)f)('n(li\ A, (), 



lns(H'tin<j: these 



vahies in the expression for r^// shows that b\, is iruiepemh'nt 
of forward s[)eed. 

Tlie variatio!! of control force witli iionnai acceleralion in 
a steady pulUu}), with no chantre assumed in forward speed 
(see n^erence S), can he foinid from (M(uations hy e(piatiMj: 
to zei'o all terms containinir I) except 1)0. This proce<lure 
implies (hat the normal acceleration is due entirelv to 
<'urvature of the di^dit path DO. The e<|uations heeome, for 
an af)|)lied hinjre moment, 



-^0 



If elleots of slipstream on the tail are ne^rlected. 



from which 



If the normal acceleration is nq. 



De 



and 



dn 



These formulas for and are equivalent to equation (i) 
of reference 0 and ecpiations (27) and (28) of reference 10. 

The formulas indicate that the stick-force ^rradients F„ 
and F„ are dependent on most of the aforemention(»d airplane 
and elevator parameters. Fijrures i to 5 show the variation 
of these stick-force <rratlients witii the parameters C,,^^ , {\, 
(\n^, //. and /i. The <i:radients are independent of spcu'd. 
altfiouirh only within the limits of tin* assumptions nuide in 
the analysis, namely. lU'i^lect of power and of eompressihiiity 
ellVcts. The ^n*adient Fu ciui l)e usi'd to ^rot the stick force 
for only a small chaiiire in forward speed f>ecause the stick 
force is not directly proportiomd to the change in speed. 
The stick force in a steady pull-up /\, however, is [)ropor- 
tiomd to th»» normal acceleration provided tlie control de- 
flection is not so ^'reat that the Imsic assumption of lim^arity 
is violate<l. 

The Iin<' () is the houndary for true static stahility — 
(hat is, Fu ^0 is the condition for zero variation in stif-k 
force with forward speed in steady (liirht. This con<lition 
is the same as that ohtained hy s(»ttinfj^ F—i), where F is 
(he constant term of the si.xth-order stability e(iuatiou 
obtained from e(|uations (1). On subsequent fitjures it is 



called the divergence boundary. The line F„ = 0 is the 
boimdary for apparent static (or maneuv(M-injr) stability and 
IS the (condition for zero variation in stif-k force in a steady 
pull-up. This condition for i^„ = <) is obtained by setting 
F={) in the approximate stability equation (equations (3)), 
which is for three degrees of freedom (a, De, and 5,). On 
the unstable side of F^^i), a slow divergence occurs that is 
noticed by the pilot as an unstable variation of stick force 
with forward si)eed. The stick foi (h> due to normal accelera- 
' tion in a pidl-up is stable, however, unless the conditions are 
; such that the airplane is operating on the unstable side of 

Figun^s 1 to 5 indicate that the parametei-s have the same 
(Effect on and F„ exet^pt that the altitude alFects only F,, 
They show that tlie stick-force gradients on an airplane of 
given tail size and center-of-gravity position may l)e in- 
<'reased by making the floating tendency (\,^^ more positive 
or by m.iss unbalancing the eh^vator control system to 
depress the elevator (make it tailheavy). The effect of the 
resloririg tendtwicy (\^ „a the stick-forc.' gradients depends 
on the relative position of the center of gravity and the 
aerodynamic center. If the c(>nter of gravity is ahead of 
th(» aerodymimic center (airplane stable with controls fixed), 
increasing the nuignitude of C,. incn»ases the stick-force 
gradients. If the center of gravity is behind the acMody- 
namic e(Miter, this effect on F, is nn-ersed; the effect on F„ 
IS not reversed, however, until the center of gnivity is well 
behind the H<M()dynamic center (in this case, about O.O.V at 
sea h'vel and 0.02c at 30,000 feet). If r,^==(), the stick 
forces are independent of the position of the airplane center 
of gravity. 
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(U'poudinir on 
lino in (iirmv T) 



Mltitiulr will (MtluM- ii\rmisi« or iiiMtTiiso 
tho hiiiirc-inoincut panunoliMS. Thr solul 
isllir locus of vnlurs of r V,, aiuU V, tof ^viu(■ll 
F is iiuUiMMulcM.l of MltiliHlc. For points to th. Vft of this 
li,;, F„ <Urmisos with nltitu<h«; for i)oints to th. n-ht o 
this lino, F„ incmisrs witii altitude. This lino is <iotrr.nimMl 
by tho rohition 



' "5 



which, for the case of lij^nirc 



f). become! 

--=l.r,OfA, 



\nothef iuct.l.o.l of inereiisiii!: the sli. k-foire irnul.e.it in 
level llisrht t\ .'onsists in applvin- a eon.Umt lunsre mo.ne,, 
to tl.e ei..vator by nu-ans of a si)rins,' or bun^'eo. The .'l e.'t 
of the sprinsr »^^ the ^'n.dient t\ is .luc to the .hTivative 
(\ whi<-l. <lepoiuls in the- same way on th.- eonstant hinsro 
.nomeut. whether it is .•ausecl by a weij:lit or by a spmig A 
,,„„,,oe, which tends to depress the e evator, w,ll thereiore 
i„ere.ise the stiek-foree gradient in level flight K- Hit 
elFect of the bungee on the stick-force gradient ,n accelerated 
flight Fn will be /em because its action depends solely on 
changes in forvvar<l sped. Its elVect on the short-per.od 
oscillations will be zero for the same reason. 

DYNAMIC STABILITY 

NO t-RICTION IN CONTROI, SYSTKM 

Trie stabilitv of the short-perioil oscillations without 
friction is sho^n in figures (i to II, which also show the 
bouiularies for true static stability (divergence '""'"•l;''-"- ; 
Ki.riuT (i is an e.vample of a more n.>arly complete presentatto. 
of^he stabilitv .lata than subse.iucnt figures b.-<ause ,( 
shows the variation of damping and perio.l of osc.Ua ion 
with the hinge-moment patameters and f for certain 
(ixe.l values of the ..ther para.net.MS. The damping, winch 
is proportional to |, increases with the niagnitu.le of ( 
The p.-iio<l, iiropoitional to ^, decreases as (\,, increases. 
\noth.-r wav of presenting this ail.Utional stability data is 
;iu.wn in ligure 7, which gives the number of cy.;ies the oscil- 
lation performs before it damps to half a.nphtude It is 
elear from ligure 7 that the osc^illation is very wd .lamped 
„„l,.ss the iTstoring tendency is .•lose to /.ero. In tins parti- 
,„,„, ,„se, onlv one oscillatory mo<le exists. Inasmu.'h as 
there are onlv tlwoe roots in this .'ase (bwause an.l <,-( ). 
,1,, other root is always ••.■al and is of no parti.'ular signili- 
eaiice in dvnamic stability. In cases in whi.-h an a.hlitional 
oscillatorv mo.le exists, it is always stable. 

The elle.t of center-of-gravity position on the =^1" "l't>'" 
the short-p.Mio<l .)s<'illations is shown m figure 8. Tlie slult 
i„ the dvnami<-stability boun.laiy, for the .•omparativo y 
hiP'e shift in center of gravity shown, is practicallv negligible. 



\r,„iv of the sul)s,-<|nenl, liguiTs, n, win.-h /..-ro static stability 
is assume.! lo fa.-ilitate .•ompntalion, th.-ielore are valul lor 
airplanes iiaving a margin of static stabiUty. 

Til.' .'(lV< t of moment .)f inertia .if the .'l.-vator <onlr.)l 
system on the .Ivnamic stability is shown in ligure i». whi.-h 
.Mv,>s tvpical vahies .if the numu tit of inertia. Tlie <''l.-<;t 's 
slightly (lestal)ili/.ing especially for high values of C,„^. I lu- 
d.-stabili/,ing .Hect of th.- inom.-nl of inertia of the .-levator 
is .M-eater than that of th.> monient of inertia of the contro 
stick n.-caus.- the ac.'uia.'V gauu^.l by in.-lu.lnig moment 
of inertia is small ompare.l with th.' saving m labor game, 
bv neglecting it. moment of inertia of tli.' .-levator .•ontrol 
svst."in was set .'.lual to zero in the subse.iuent calculation. 
As a ivsult, the stability equation Ix-comes a cubic and suli- 
sequeut ligures are somewhat unconservative. 

The eir.>et of .lensity iiaraineter n on the dynamic statiUity 
is shown in ligur.- 10. " Increase of m has a slight destabilizing 

'"'vs 'has be.-n shown, mass uniialan.T of the elevator contml 
svstom improv,.s the static stai,ility dig. 4). The ellcct on 
.iynamic stability is unfavorabl,-. howev,-r. as shown in 
ligure 11. The value of mass unbalance shown corresponds 
to a bobweight that is located at the airplane center of 
..ravitv and ie<(uin-s a balancing pull of ;57 pounds on the 
control stick of a i,u.-suit airplane at sea level. Increasing 
oscillations o.cur if the aerodynamic balance is too high 
(low magnitudes of C,,) , especially for n.-gative values of C ... 

The eirect of the location of the bobwcight is shown in 
ligur.> 12. Kach curve represents a .liffcrcnt arrangement 
of bobweights and all arrangements give tho same stick 
force Tho solid lino corrcspon.ls to a weight at tli.^ air- 
plane center of gravity (for practi<-al purposes, at the control 
sti. k) The short-da.sh line corresponds lo a weight at tlie 
elevator. The long-.lash lino correspomls to two weights- 
one at tho elevator, which ton.ls to make it noselieavy; the 
other at the .-ontrol sti.-k, which gives the elevator a sulh- 
. ieiitlv powerful tailhoavy momem that the resultant stick 
f„r(-.- is » lie same as wit li the single weight . In the particu- 
lar case repiesent<-.l, the noselieavy moment duo to the 
weight at tho <-l.-vator is etpial to the tailliea%-y moment due 
lo both weights. Moving tho single weight from the con- 
irol stick to the ,-levator has a large ,l(-stabilizing eirecl. 
Overbalancing the elevator while tl..-. stick !ovca is kept 
.-onstant has a .-onsi.lerablc stabilizing eirect. This metho. 
of prev-nting instability has the .lisadvantage. however. .)f 
in<-irasing the total amount of unbalancing weight requir.Hl. 
In the ,-ase shown. Ihe weight is incr.-ase.l to three tiin<;s Us 
original si/.e. Another dismlvantage is the roarwar. shift ... 
..enter of gravitv of the whole airplane .hie to additional 
w.-ight at the elevator. (See airplane paramete.-s given in 
"Method of A.ialvsis.") The dcstabili/.ing elfecl of the i.i- 
e,-eased monient of inertia asso.-iated with the ml.lod weights 
was fouml to b.- v.-ry small, especially for negative lloatmg 
tendency. 
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KI-FK( T OK VISCOrS FRICTION IN CONTROL SYSTEM 

In (ho pivcc^liiiii: seel ion. a constjuil vnliio of IIk^ olovntor- 
(ianipinii: parainolor (\,^^^ was assuinpd. This valno was diic 
only to aorodyntiniic (hinipinir. Tho cflVcts of viscous fric- 
tion in th(» (»I('vator control system luv oi>t}n.nc(l by consi(h»r- 
iiiir as an achlitional variable. This vsirijiblc <'an l)c 

introduced. Jis in the preeedin<r section, by showini,^ a. seri(^s 
of boutuhuies in the ^V,^^^ plane foi' various vahies of 
C Tlie iretu'ral nature of (he (^H'eet of friction is shown 
firs(, however, l)y piHsentin*: b()un(hiri(»s in tlie (ji^^^/.j,^ 
phuie wilh G,,^ constant and some otlier paranu^ter varied. 
This method of presentinir stabihty boun(hiries mak(^s it 
easier (o show (he etfiM'ts of otluM- [)arameters sucli as air- 
phme center-of-trravily position and (h'lisity wiieii friction 
is iiitro(ht(*e(i. 

The etfeet of viscous fri(*tk>ii on tlie <lynamie stability, 
for various conditions, is shown in fiijures \:\ and 14 for 
M=12.r) and ^ = 37.5, respectively. Fi<]^res i;Ua) and i4(a) 
refer to the mass-balanced elevator control system; figures 
K^b) arid 14(b) r(>fer to the tailheavy elevator control 
system coiisuiered in tlie precedintj; section. It is shown tliat, 
if the airplane center of jrravity is ahea<l of a eortaiu point/ 
tiie instability caused by the unbalanced elevator can bo 
removed by a<lding viscous friction to the control system. 
This critical center-of-gravity position is behind the aero- 
dynamic; center, an<l its <[istance from the aerodynamic 
c(*nter decreases as the density parameter m increases, 
Wh(»n the center of <j:ravity is behind tliis critical position, 
viscous friction has a (h^stabilizing effect. These opposite 
efFects of viscous friction are shown in the C^^^ Cf,^ plane in 
fifjjures If) and 10. AVhen the center of gravity is slifjhtly 
ahead of (his cri(ical position, the effect of viscous friction 
depends on its magnitude and also on the value of C^^ The 
additioii of a small amount of viscous friction is destabi- 
lizing but larger amounts are stabilizing. If the aerodynamic 
balance is sufhciently high and the viscous fri{'.tion 

lies in a certain range, increashig oscillations will o<'(;ur. 
Jn figure 14(1)), for example, if — O.OIc and G^^— 0.05, 

the oscillations will be unstable when the elevator-damping 
parameter is in the range from —2.5 to —70. If is 
mon» negative than — O.OSf), no amoimt of elevator dump- 
ing can cause increasing oscillations. As the center of 
gravity inovcs forward, tlie destabilizing effect of elevator 
damping b(»comcs less and finally disappears. 

Tile effect of tlx^ density parameter ^ e»n be seen by com- 
I)aring figures 13 and 14. The critical center-of-gravity 
position approach(»s the aerodynamics center as m increases. 

' Siru'c I Ins n [M)ri \Mis wi if ten, llii.s|)(»iiit liiis Immmi fouinl to Iw whrro 0,.s(>in(»iinu'f!('aiU»<l 
(lu> St ifk-ilxr<l tiiiuu'tiMT (loirit. 



When M - 12.5, elevator damping always has a stabilizing 
effect provided is f)()sitive. When /i = :^7.5, elevator 

damping may be destabilizing ov(m- a small range of C,,^^^ 
and (\,^ even when .r,,.,. is j)ositive (0.05c). 

When the center of gravity is slightly alw^ad of the afore- 
mentioned critical position (w^iich is behind the aerodynamic 
cent(M), the conditions under which elevator damping may 
cause dynamic; instability may be advantageously repre- 
sented in the C\^C,^ i)lane. If a series of stability boundaries 
are drawn in that i)lane for various values of elevator dam])- 
ing, they will all be confined to a region bounded by a line 
that will be called the boundary of com])lete (hunping. An 
illustration of two methods of constructing this boundary is 
given in figure 17. If a series of !)0undaries in the Ch^Qj. 
I) lane are draw^n for various values of the damping, the com- 
mon tangent of all these curves is the boundary for complete 
damping. This boundary can also b(* drawn by plotting the 
minimum valuers of C^^ obtained from ])lots of the type shown 
in figrnvs and 14 aganist coiirsponding values of C^^ , 
The region in the C^^Cn^ plane between the boundaries for 
complete damping and increasing oscillations is the region 
where the addition of viscous friction to the elevator c:ontrol 
system may cause dynamic instability. 

That a boundary for complete damping cannot be shown 
for M= 12.5 if the airplane is statically neutral or stable (x«.c. 
is zero or ])ositive) may be seen from figure \\\. It is possible 
however, to show a boundary for complete* damping under 
these conditions for /i = ;i7.5. Figure 18 shows these bound- 
aries for a:«.tf.=0 and for the critical value arfl.c.= — 0.017c, 
for both a mass-balanced elevator and a mass-unbalanced 
elevator. The boundaries for increasing o.sc illations and 
divergence are also shown. For the case of the mass-balanced 
elevator (A = 0), the boundary for complete damping is a 
straight line through the origin and thcM-eforc^ corresponds 
to a fixed ratio of the floating and n»storing tendencies, 
or floating ratio. Elevator mass unbalance decreases the 
region of com])lete damping. 

EFFECT OF SOLID FRICTION IN ELEVATOR CONTROL SYSTEM 

The boundary for com])lete dam])itig has an important 

bearing on the effect of solid friction on dynamic stability. 
In order to calculate this effect, the solid friction is re])laced 
by an ecpiivahMit viscous friction that would dissi]}ate energy 
at the same rate. This condition gives an eriuivalent 



for a sinusoidal motion of the elevator with amijlitude 6 and 
angular frcHpiency 77. 
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KiOL'RK 17.— Method of c<)nstructinK boundary for coinplcto dampittR. 



It rail he shown that if vis<'Ous friction is <[rstal>iiizinir, as 
in fiij;uros Hi (o 18, solid friction may lead to strndy oscilhi- 
lions havinj? an amplitudo ])roportionul to the aniouut of 
friction. Suppose an os<'iihitton is start('<i with ampHtude 
and froquoncy wliirli result in an e(|uival(Mit C^^^^ that would 
cause inereasin«r oscillations. Let this conditioti he r(*])re- 
sented hy point 2 in fii;:urc 19. The amplitude of the oscil- 
lations would then increase until the e(juival(»nt C\^^ de- 
creased to the value that would result in neutrally damped 
oscillations. repres(>nted h}^ point W in fi<^ure H). If the 
initial amplitude is so low that the ('(pavalent viscous fric'- 
tion is in tin* slahle re<^ion, as shown hy point 1, the oscilla- 
tions will die out completely. If the initial aniplitud(» is so 
hiirh that the osi'illations are stahle, represented hy point 4, 
the amplitude will decrease until it reaches a constant value, 
when the e(pnvalent L\^^ is apun re])resente(i hy ])oint M. 
The value of at point :{ then determines the amplitude 
of the steady oscillations. By sol vino- formula (4) for 1, the 
am])litude of the steady oscillation is ohtained as 

wlu're 7? and Chj^^ are the values at ])oint W. This fornnda 



shows that the amplitude is ])roi)ort iomd to the amount of 
solid friction. 

The fore^roinir analysis shows that the repon in the Cn^Ch^ 

j)lane hetwecMi the houiulary for itK'reasinir oscillations and 
tlu^ houtulary for complete (lamping is the re<rion where 
steady o.scillations nuiy occur hecause of the induenee of 
solid friction in the control system. All the remarks in the 
])recedinij: section con<'ernin<i: the houndary for complete 
dampin<^ with viscous friction <'onse(|uently apply to the 
houndary for steady oscillations with solid friction, inasmuch 
as these two houndaries are the sanje. within the limits of 
the assumptions involved in the use of the concept of equi- 
valent vis<'ous friction. Steady oscillations due to solid 
friction will not occm* on a statically Uiuitral or stahle air- 
])lan(\ for instanc(\ unless m \^ vc^-y larp» (corresponds to a 
hi«rh altitude). Kv(mi in that case, the ])ossibility of steady 
oscillations <\xists only for a comparatively small ranire of 
float iuir ratios. If tlu^ airplane is statically unstable hy a 
sufileient amount, however, steady oscillation may <»xist 
ov(M' the entire ranire of floatinjj: ratio. 

Some calculations of the amplitude of the steady oscilla- 
tions, e.\press<'d in terms of normal acceleration per imit 
frictional force as f(4t at the control sti<'lv, were made bv the 
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f^estorin<j tendency, 

(a) A -10. 

(b) A-0. 

Figure W.-Boundaries for complete dampinK, increasing oscillations, and divergence 



method of appendix The results are presented in fi<?uro 20, 
which shows lines of constant amplitude in (he C„^ C^^ 
plane for an airplane with the center of trravity at the critical 
position ref(Tred to in the preceding: section. Steady oscilla- 
tions will therefor occur throughout the entire rejrion whore 
stability exists in the absence of friction. The amplitude 
is smallest for hiirh values of Cn^^ combined with hi^di 
values of C^y 

CONCLUDING REMARKS 

The stick-free static stability of a conventional airplane 
may be improved by making the elevator floatinjr tendency 
more positive or by mass-unbalancin<r the elevator control 
system to make the elevator tailheavy. lucreasinir the 
restorinsj: tendency also has a favorabh* ellV'ct provided the 
airplane is stable with stick fixed. If the restorinir tendency 
is zero, the stick-free static stability is independent of air- 
plane center-of-irravity positiotK 

The dynamic stability with frictionless controls <le|)en(ls 
c liiefly on the restorin^r tendency (\ and on the mass balance 
of the elevator control system. If the <»levator control 
system is nuiss unbalanced (elevator made tailheavy) by an 
offset weijrht at the control stick and if the restoring,' tendency 
is too low, inereasinir short-period oscillations nuiy restdt. 
This condition can be remedied by the use of two additional 
weitrhts— one at the elevator makin.ir it noseheavy. the other 
at the control stick making the elevator sufficiently tailheavy i 



t hat the combined effect gives the elevator the desired amount 
of tailheaviness. 

The addition of viscous friction to the control system will 
prevent dynamic instability provided the airplane center of 
gravity is forward of a critical position which is behind the 
aerodynamic center and approaches it as the value of the 
density parameter ju increases. If the airplane center of 
gravity is Ix^hind this critical position, viscous friction will 
have a destabilizing ciYvct, no nuitter what the hinge- 
moment paramet<>rs are. If the eenter of gravity is slightly 
ahead of the critical position, viscous friction may be de- 
stabilizing for a limited range of values of viscous friction and 
the hinge-moment parameters. A low rc^storing tendency 
and a high positive floating tendency will tend to cause this 
dynamic instability. When ^ is very large (high altitude), 
this condition of steady oscillations can occur even if the 
(•enter of gravity is ahead of the aerodynamic center. 

The presence of solid friction in the conlrol system may 
cause short-i)eriod steady oscillations under the conditions 
for which viscous friction is destabilizing. The amplitude of 
the oscillations is proportiomU to the amount of friction 
present. 



LaNGLKY .NfKMOKIAL AkKOXA T'TICA L La nOKATOIi V, 

N'ational Advisory Committf:]-: for Akroxautics, 
Laxgley Fikld, Va,, Dccembfr ^3, 1043, 
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APPENDIX A 



EVALUATION OF STABILITY DERIVATIVES 

Derivative . -T\h) total liin^re m(»inont acliiiir on the wIum'O 
t'lovator may Uv ox[)rcsscMl as 



At trim //=(). tlioreforo. 

lly= — P V T ' 

If //o= — 4 — » 

Using ;^ pF2CL5„=m.7 gives 

^ _ ^cg pCLSt, _ _ ]}:0l_ 
^'*"~* 2mf7 2.4«.ju 



Derivative r,„^.-Tho parameter (?..^ may l)o obtained 
from wiiui-tunncl measurements or, if the position of the 
aerodynamic center of the complete uirphine is known, may 
1)0 calculated by the formula 

Derivatives (\^^ and ^^^^.-Thc derivatives T.^^ and 
C,^ , arise because of the la^ between the chanire m an^de 
of "attack at the wing and tlie corresponding downwash at 
the tail. It is assumed that the downwash at any instant t 

depends on the angle of attack at the instant t - yy the dilTer- 
ence being the time required for the air to move from the 
wing to the tail. If a=/(0, tins relation may i>e expressed 
as 



Now, 



Ife 



nee. 



fit^ AO -y«) - (t) r (0 - 



or, because «==^ ^« and a= -q^ D'a^ 

and, because a/=a— «, 

and 

The part of the pitching-moment coefficient contributed 
by the tail is 



The lag elFectively introcluces derivatives Cm^^j ^^i^a^ . . . . 
The first two of these derivatives are 



and 
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Derivatives C^/ ^-^"p.' "'^^'^'^ (U'rivntivos 

r,,^,^^, luid C\j^>^ may l)o ()l)laiiica from 



wliicli *rivcs 



Derivative -Tho i)iLchintx nioinent (lm> to pitchiiis: 

is mmli' up of parts duo to ])ro!>('llor, ^^inJr, fusolajjo, and 
horizontal tail. Tlio conlrii)uLio!i of tlic tail is by far tlir 
largest and can oasily lu* caloulaUHl. 

If tho airplane is rotatinj: willi aujrular vrlocity 0, tlio 

inoreaso in angle of attack at tlie tail is ^ ''^^'''^^^^ 
in an increased lift on the tail given (in coeifieient form) i)y 

The resultant pitehinj^-moinent coodicient is 
tmd cxpressini? 0 as ^ 06 and " as //, gives 



Derivative -The derivative r,„3 may he measured 



<hreeLlv in a wind tunnel or mav l)e eomputed from wind- 
tunnel (lata on the vahie of (\,, for tlu^ horizontal tail by 

means of the formula 



r 



Derivative 6',.,,,,. -Tiie derivative C',,.,,, nniy he eomputed 
from 

(oS. ''''' '''''' ^ 

<»f reference IK vvliieh is based on thin-winjr potentud-llow 
tiu'ory. 

Derivative C,,^^^. The derivative^ (\,,^, is given by 

In tlie absence of viscous friction in the elevator control 
system, the value of C,,,^ may be comi)uted from 

whore (^^Js), (dm),, '""^ ^ 
of reference 1 1 . 

If a dashpot, whidi has a (himping constant of K pounds 
per foot per second and moves a distance of Q feel per 
radian of elevator dellection, is inserted in tlie control 
svstem. 



(A2) 



The contribution of the wing depends on the assumed 
axis of rotation (center of gravity) but a fair average vahie , 
will b(» obtained 1)V assuming tliat the center of gravity is i ^^i^^ t^tal value of T,,,, is the sum of e(|uattons (Al) and 
at the wing ciuarter-chord point. This assumption gives a (X2) 
value 



The U)tal pitching-moment coefficient duo to pitching 
therefore is 



Derivatives T,^ and C,,. -The derivalives T,^^ and 
can be calculated by thiu-wing-section theory but the 
results are of doubtful accuracy because of three-dimensional 
and boundary-layer etVects. It is tlierefore best to obtain 
these derivatives from wind-tunnel tests. 



(UI) 



APPENDIX B 

< ALCH,,AT.ON OV NORMAL ACCELERATION DUE TO OSCILLATING ELEVATOR 

ulual t., IHa^e^ u, j.o,ul,.nc„s.(,„«l units, ca,i I... cnlculatcl from 

/J(a—0) (. ^ 

" I r "r~r W-^- ^-^ 

•) I "'ne^^a ' — ' 

wJuMV 77 is {\w ajiiriilar fmiiK^nry of I ho olrvator. 

^«HU <Mii i)( (onv(Mto(l to i)hysicai units by tho formula 

Xorinal acceleration por (/ ^ 4/,r /)(a—0) 4 
Stick friction in pounds' ' pS^c^g " ^ ^^776', 

wiu>re C.,, . the value of elevator <lampin, requiredjor the eon<iition of neutral dynanl stability. 
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2, ccl-a X, sl^ctb una fr« bottcM. The ^sBt-g te» on th. left- 
MM side of the equation ehould be ^s- 
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Force 
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^Mdment about azw 



De8i^tio£ 
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PfiOPELLEB SYMBOLS 




1 hp«76.04 kffVs=*^50 ft-lb/sec 
1 metric horsepower^ 0.9863 hp 
1 mph=0.4470 mps 'Z^":--' ..i^r- 
I mp9=i.23«9 mph- : ; ^ ' ^ , 



Mb=0.4536kg 
1 kg«:2.2046 lb 
1 nri*= 1,609, 36 ia:^5;280 ft 
tni«3;2808ft V ^ - 



' ..... ^X'"^^ ■'' ■. 



